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ABSTRACT 



The major environmental Ultra Low Frequency (ULF), magnetic 
noise sources observed on the airborne AN/ASQ-8KV) optically 
pumped magnetometer are electromagnetic micropulsatrons , geo- 
logical structures and ocean wave noise. 

Experimental determinations of the amount of magnetic noise 
produced by ocean surface gravity waves were made. The author' 
current research designed to remove the wave noise and increase 
the probability of Magnetic A.nomaly Detection (MAD) is summariz 
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I. INTRODUCTION 



Magnetometry has been used as a classical geophysical pro- 
specting tool for years (Dobrin, 1952 ). The principle^ o.f 
magnetic anamoly analysis has also been applied to Anti- 
Submarine Warfare (ASW) with a limited amount of success for 
submarine localization. 

The sensitivity of magnetometers has been continually 

improving so that it is now possible to observe, from an air- 

_ 2 

craft, fluctuations of 10 nano Tesla/hertz (nT/Hz), imposed 
on the earth's constant field of approximately 50,000 nT . 

Land stationed superconducting magnetic Sihesors have demon- 
strated a sensitivity of 10~ nT/Hz. The advent of airborne 
superconducting magnetometers would greatly increase detection 
ranges 5 however, the extent of improvement depends in detail 
on the amplitude and frequency of the environmental noise which 
is above the noise floor of the instrument. 

Airborne Magnetic Anomaly Detection (MAD) for ASW is done 
in the .04 to 2 Hz frequency band. Ocean surface gravity waves, 
in the .03 to 10 Hz frequency band (Silvester, 1974) are con- 
tributors to the magnetic noise spectrum. The other two major 
environmental noise sources are electromagnetic pulsations and 
geological noise . 

It is generally assumed that noise due to ocean gravity 
waves is not "loud" enough to be of any consequence. However, 
experimentally it was observed that noise levels induced by 
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typical sea condition in the mid latitude oceans during the 
late fall to early spring are at least one order magnitude 
greater than the noise floor of airborne magnetometers. The 
noise greatly reduces the effectiveness of Magnetic Anomaly 
recognition . 
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II. THEORY 



The problem of determining the electric current induced 
by deep water gravity waves at the ocean surface and the resul- 
tant periodic fluctuations in the magnetic field has been de- 
scribed by Weaver [1965]. The magnetic field fluctuations 
above and below the ocean were obtained by solving Maxwell's 
Equations for the magnetic flux density vector in air and water, 
and applying electromagnetic boundary conditions at the sea 
surface . 

The movement of the sea surface through the earth's mag- 
netic field induces electric eddy currents in the sea water. 

In deep water, assuming linear ocean wave theory, the water 
particle motion is periodic and circular, as is the induced 
magnetic field. This process becomes more complicated as one 
decomposes the sea surface spectrum into waves of various 
frequencies which may be propagating in various directions 
with respect to the magnetic field vector. Further complica- 
tions arises as higher order wave theories are applied to de- 
scribe the water particle motion and the waves propagate into 
shallow water. 

A. MATHEMATICAL MODEL 

The mathematical model of the induced field was developed 
by Weaver [l965], assuming linear wave theory to describe the 
water particle motion. 
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It was derived for a monochromatic ocean surface gravity 
wave, propagating in deep water, in the X direction. The 
development is summarized below. 

The coordinate system formed a right hand triad as follows : 

X-axis, was taken in the direction of wave propagation. 

Y-axis, was taken parallel to the wave crest. 

Z-axis , was taken downward into the sea, Z=0 at the 
still water level (SWL). 

It was assumed that the earth's magnetic field was constant 
in the region, and the wave crest was infinitely long. The 
model became two dimensional as all magnetic gradients in Y 
were negligible. 

The earth's vector magnetic field, 3^, was expressed as 



where i, j, k were unit vectors along the three coordinate 
axis X, y, z. I, was the magnetic dip angle in the vertical 
plane, relative to the horizontal. 9, was the inclination 
in the horizontal plane of the direction of wave propagation, 
clockwise from the magnetic meridian (see Figure 1). 

The field vectors satisfied Maxwell's equation in the form 





( 1 ) 



( 2 ) 
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Figure 1. The model as it applies to a one component sea. 



where was the electric field vector (^°^'*‘®/Meter) , B was the 

induced magnetic flux density (Tesla) , was the permeability 

of free space ( 4 ttx 10 ”'^ ^^^^^/Meter) , o was the electrical con- 
Si0iTi0n 

ductivity ( /Meter), V was the water particle velocity 

vector (^®^®^/sec), e was the permittivity ( ^ ^^^'^^/Meter ) , and 
t was time (seconds). 

The deep water particle velocity, described by linear wave 
theory, was of the form 



^ 63( Cc-f ^)i2. 




(3) 



where a was the wave amplitude (Meters), to = 27r/f was the angular 
frequency of the wave (^'^'^^^^^/sec) , m=2Tr/L was the wave number 
(Radians^jl^^g^^ ^ and t was time (seconds). 

Since V was harmonic in x and t, the magnetic and electrical 
field resulting from V has a similar structure, thus 

L (_uJC-/V»lX) 



?> ' bc-^) iL 
E ^ E(0) z 



lC^IC. -/wX) 



(4) 



(5) 



The B field and the E field are non-divergent which 



— LffA D A 



-z L m t-vC-3.) 



^ 3 ^ 



implied 

( 6 ) 

(7) 
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The differential equation in water satisfied by B was 

Ji, Ok 

obtained by assuming B<<B^ and coe<<a. Then, taking the curl 

of equation (2) and using the non-divergent portion of the 

^ 2-^ 

vector identity, Vx(VxA) = -V A + V(V-A), the following was 
obtained : 

- V^5 = G' vxf -f e_/{oVx(v:< b^) . ( 8 ) 



Further simplification of equation (8) was obtained by the 

^ k 

elimination of VxE by equation (1). From equation (3), V was 

2 

eliminated. Then if r -- iy cocf, the following equation results , 






9) 



In region Z>0 , equation (9) was solved for the Z component 



by assuming 






where Q was an undetermined coefricient. 

In air a>>a, thus equation (2) became, 

vx-S- ^ 

Ao (yt 

^ . 



( 10 ) 



The curl of elation (10), when E was eliminated by equation 



( 1) , resulted in 






( 11 ) 



Confining the sea surface to scales of motion described by 

2 2 

gravity waves implied that m >>y^ea) , thus 
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In the region Z<0 the solution for the Z component of 



equation (12) was of the form 



( 13 ) 



where R was an undetermined coefficient. 

The magnetic field vector was continuous across the air-sea 
interface, thus, the following boundary condition at Z=0, existed 



U-ATeC. 



(14 .a) 



Alt (2=o) 















(14 .b) 



Applying the above boundary values to equations O) and (13) 
yielded . . \ 






(15) 



and 



= -4*64 /»i^) - - 1 'Sic.j)) , ( 16 ) 



Equations (15) and (16) lead to 



/>n i-(^ 



( 17 ) 



and 






( 18 ) 



Thus the Z component of the induced magnetic field vector 
in air was 



■>‘1 C^iCx) ■"/^) £. ^ 

" . 2 — ^ 



9) 



In order to further simplify equation (19), — 2 "^^^ true 

m 

for gravity waves resulting in 

From equation (6) and (20) it was seen that 

ib^C^) (21) 

It was convenient to measure the earth's field with a 
total field magnetometer. Thus, 

= B^+ h,C^) Calces &+b0, S,Aj I . (22) 



When equations (20) and (21) were used to simplify equation 
(22), and only the perturbations were considered, it can be 
shown that the magnitude of the observed magnetic field at a 
specific frequency would be. 



■k/')= ^^l^i((ccsxccef^(^,^ xf) rS' 



y 

z<c 



(23) 
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The deep water wave length (L) of a surface wave is related 
to the wave period (T) by, 

9 "^^ (24) 

"TTr 

where g is the acceleration due to gravity. 

Eliminating m, from equation (23) by the use of equation (24) 
leads to, ^ 

0*1 • 

Equation (25) was used to model the magnetic field spectriom 
when the directional ocean wave spectrum was known. 

Using an airborne magnetometer to measure the magnetic 
spectra, there was a doppler type frequency shift, described as, 

- ( Zrrf^/^)coc^) 

(see Appendix C). Where f was the wave frequency observed on 

the aircraft (sec~^), f was the actual wave frequency, u was 

Ms "t 02T* 

the velocity of the aircraft ( /sec) , and was the angle 

between aircraft morion and direction of wave propagation (see 
Figure 1) . 
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III. THE EXPERIMENT 



Magnetic data was collected by an airborne magnetometer 
flown over the sea surface in order to verify the theory pre- 
viously presented and to further understand the observed environ- 
j ment in the 0 to 2Hz band. The AM/ASQ-81 magnetometer, aboard 

I 

the P-3C Orion aircraft, is capable of measuring the signal 
produced by ocean gravity waves which are frequency shifted 
due to relative motion between the wave and the aircraft flying 
at 100 meters altitude. The analog voltage time series output 
of the magnetometer was recorded on an HP3960, FM tape recorder 
for post flight analysis. 

The fluctuation of the sea surface below the aircraft was 
determined by two techniques. The analog voltage time series 
output of the APN-194 radar altimeter, also aboard the aircraft, 
was recorded. 

Independent of the aircraft data the ocean surface wave 
climate was determined from analysis of the Fleet Numerical 
Weather Central (FNWC), spectral Ocean Wave Model (SOWM) and 
measured wave spectrum from National Oceanic and Atmospheric 
Administration (NOAA) , Environmental Buoy, NR-6003 (E3-17), 
located at 52.2N, 155. 8W from the airborne data. 

o o 

The airborne data was collected in vicinity of 52 N, 156 W, 

I 

around the buoys' reporting period, OOZ, 10 Feb 78. 

A. AIRBORNE DATA COLLECTION 

The AN/ASQ-81 employs the atomic properties of optically 
pumped, metastable helium atoms to detect variations in total 
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Figure 2. Region of data collection. (X) FNWC Grid Points, 

(A)NOAA Buoy, ( ) First Leg of Aircraft Data, 

(_ ) Second Leg of Aircraft Data. 
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field magnetic amplitudes. The output of the sensing element 
is converted to an analog voltage which is bandpassed and scaled 
to display field variations (Order 10~^^T) as a function of time. 
Instrument specifications are described in Appendix A. 

Perpendicular legs, 270*^ True and 000*^ True, were flown at 
approximately 60 meters altitude and 180 knots (91.44 M/s) for 
30 minutes each. (See Figure 2) The magnetic data was collected 
in the .04 Hz to 2 Hz band, thus, eliminating the earth's con- 
stant field and any micropulsations above 2 Hz. Aircraft motions 
were compensated for by the magnetometer using baroalt imetry , 
selective filtering and accelerometers . No variations in the 
magnetic field due to the aircraft were expected due to the 
above real time processing and the straight and level flight 
paths flown . 

The instantaneous aircraft altitudes were recorded simul- 
taneously with the magnetic data from the APN-194 radar altimeter, 
which responded linearly to apparent altitude changes due to 
long period surface gravity waves . The wave frequency was 
shifted due to aircraft motion (see equation (26)). Beam 
characteristics, 6 0*^ cone, were such that at 6 0 meters the 
diameter of the radar spot on the surface was approximately 70 
meters. In order to resolve wave amplitudes from 60 meter 
altitude, the wave period must be greater than 10 seconds. 

The frequency and amplitude errors attributable to the size 
of the radar spot were considered in the analysis. 

The AN/ASQ-81 and APN-194 signal pick-off points are 
described in Appendix A. 
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B. SPECTRAL WAVE DATA COLLECTION 

The experiment was conducted around synoptic reporting 
period OOZ, 10 Feb 1978 in vicinity of NOAA Buoy, EB-17 , so 
the spectral information computed by the Buoy for the OOZ time 
period was useful in interpreting the aircraft data. 

In addition to Buoy data the SOWM output was used in order 
to infer directionality and improve the spectral wave estimate 
in the region of the flight . 

The Buoy is designed to report meterological and oceano- 
graphic conditions at the sea surface every 3 hours . The wave 
conditions are measured for 20 minutes by acceleromerers and 
reported for the Buoy location as frequency and amplitude 
(Steele and Johnson, 1977). The Buoy data of OOZ, 10 Feb 1973 
is reported in Table I . 



Frequency (Hz) 


Period 


( sec ) 


Amplitude 


.06 


16.7 






. 3322 


.07 


14.3 






.7574 


.08 


12.5 






1.6064 


.09 


11.1 






1.9604 


.1 


10.0 






1.9312 


1 — 1 
1 — 1 


9.1 






1.5419 


.12 


8 . 3 






1.0116 


.13 


7.7 






.9033 


.14 


7.1 






.7917 


.15 


6.7 






.6543 


.16 


6.3 






.4793 


.17 


5.9 






.4511 


.18 


5.6 






.4212 




Table I 








NOAA Data 


Buoy OOZ, 


10 Feb 


1978 
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The FNWC, SOWM uses the analyzed surface pressure field 
to determine wind velocities and then builds, radiates and dis- 
sipates waves as prescribed by wave theory (FNWC Tech Mote Mo. 
75-3). Model output is energy as a function of frequency 
and direction at each grid point. The wave energy (E) at a 
frequency can be converted to wave amplitude (a) by a= E-Af. 

SOWM frid coverage of the data area is shown in Figure 2. 
The model output for the closest grid point is summarized in 



Table II. 

DIRECTION FROM 283° 313° 343° 013° 

(°T) 



PERIOD (SEC) 




AMPLITUDES 


(METERS) 




16.4 


. 347 


0 


.578 


. 347 


14.9 3 


.578 


.408 


1.002 


.578 


13.89 


. 708 


1.293 


1.082 


. 347 


12.82 


1.636 


0 


1 .686 


1 .686 


12.04 


1 .227 


1.686 


1.157 


.959 


10.99 


1.417 


1.584 


. 914 


1.192 


9.71 


1.293 


1.493 


. 944 


1.435 


8.5 


.944 


1 .356 


1.227 


1 .271 


7.5 


.851 


.973 


.987 


.914 


6.5 


.409 


.457 


.354 


.157 



Table II 

FNWC, SOWM Data, 03Z 10Feb78. Grid point located at 
52.357N, 154.099W. Amplitudes (meters) 
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IV. DATA ANALYSIS 



A. TIME SERIES ANALYSIS 

The general form of the magnetic and altitude time series 
is shown in Figure 3 , which is a sample of each time series . 

The time series were digitized at .1 second intervals, 
resulting in a Nyquist frequency of 5 Hz. Any spectral power 
above 5 Hz would be folded back into the spectral estimate 
(alliasing) . The magnetic data, due to the low pass, 2 Hz, 
characteristics of the instrument, alliased only instrument 
noise which was at a low level. It was assumed that no sig- 
nificant alliasing of the altitude data occurred due to the 
altimeters 0.1 second RC constant which results in a 1.59 Hz 
cutoff frequency ( = ^/RC). 

Prior to computing power spectral estimates , the data was 
calibrated to insure constancy of unit magnitudes. The magnetic 
data voltages were linearly calibrated to nT and the altimetry 
voltages were linearly calibrated to meters after the flight. 

The data was prepared for the Fast Fourier Transform (FFT), 
(Bendat and Piersol, 1971) by forming 16 ensembles, 102.4 seconds 
each. Each ensemble was then linearly detrended. The detrending 
removed significant variance from the frequencies below .001 Hz. 
The cosine function was used to taper 10% of the time series 
by multiplying the first and last 51 points of the time series 
by the cosine. Points 1 and 1024 become zero and points 52 
thru 973 are as digitized. The tapering procedure is done in 
order to prevent non-real low frequency power leakage, due to 
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Figure 




Sample of magnetic time series (top) and 
altitude time series (bottom) . Horizontal 
scale is Imm/sec. Vertical scale is Imm/volt. 
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ensemble averaging, into the spectrum. The variance removed 
in the time series taper is restored to the power spectra after 
the FFT (Bendat and Piersol, 1971). 

The FFT computes the real and imaginary components of the 
Fourier Transform. It was then possible to compute the sample 
power spectra of each signal by ensemble averaging. The normal- 
ized cross power spectra, defined as coherence and the phase 
of the cross power spectrum between signals were also computed. 

Each of the two flight courses were treated as separate 

data sets. Note in Figure 4 the power spectra of the alti.meter 
2 2 

(M /Hz) and magnetometer (nT /Hz) and the coherence and phase 
as a function of frequency. The numerical results are summarized 
in Appendix D. 

B. SPECTRAL WAVE ANALYSIS 

The analysis of the FNWC , SOWN directional wave spectra 
and EB-17 ocean wave amplitude spectra in order to describe 
the magnetic field variations observed aboard the aircraft, 
required fitting the large scale information to the microscale . 
This change of scale was feasible due to spatial uniformity 
of the ocean wave spectra over the area and frequency range 
investigated . 

The ocean wave climate in an area can be described by that 
at a point over a period of time . EB-17 reported the waves 

during a 20 minute period to have the spectral shape listed in 
Table I. The FNWC ocean wave amplitude and directional spectra 
had a similar shape listed in Table II. 
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The directionality of wave components at EB-17 was deter- 
mined by group velocity analysis (Silvester, 1974) which pro- 
pagated the FNWC, SOWM data at nearby grid points to the buoy 
position. The greatest weight was on the grid point closest 
to the position of the buoy. The grid point was located at 
52.357N, 154.099W, 60 miles from the buoy. The analyzed buoy 
directional amplitude spectra is reported in Table III. 



PERIOD (SEC) 


AMPLITUDE (METERS) 


DIRECTION FROM 


16 . 7 


. 8322 


340° 


14 . 3 


.7574 


340° 


12.5 


1.6064 


315° 


1 — 1 

1 — 1 
rH 


1.9604 


315° 


10.0 


1.9312 


320° 


9.1 


1.5419 


320° 


8.3 


1.0116 


320° 


7 .7 


.9083 


320° 


7.1 


.7917 


330° 


6.7 


.6543 


330° 


6.3 


.4793 


320° 


5.9 


.4511 


320° 


5.6 


.4212 


320° 



Table III 

Analyzed directional wave data at NOAA buoy E3-17 , OOZ 10 Feb 78. 



It is important to note that the iMOAA data and the FNWC 
data have similar spectral shapes although amplitudes are not 
equal, verifying scale analysis assumptions. The buoy reports 
a spectral peak at approximately 11 seconds . FNWC predicts a 
spectral peak at approximately 12 seconds . 
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The observed ocean wave climate on the aircraft, although 
similar, was not the climate observed by EB-17 or predicted 
by SOWM. Each of the ocean wave data sources had errors asso- 
ciated with them, which would contribute to the differences 
in the data. Scale considerations, length of airborne data 
record, modeled grid point locations, buoy location and flight 
path, further contribute to the uniqueness of the separate 
data sources . 

Due to the multidirectional nature of the SOWM data and 
its proximity to the flight path, it was assumed that rhe 
closest FMWC data point (Table II) most closely resembled the 
ocean wave climate observed on the aircraft . The SOWM wave 
information was applied to equation (25) resulting in an ex- 
pected magnetic amplitude spectrum. The expected frequency 
shift due to the aircrafts motion relative to the waves was 
computed by applying equation (26). To be consistent with the 
power spectral density computed from the aircraft data the 
amplitude spectrum described above was squared and normalized 
by multiplying by the frequency interval of each frequency 
estimate. The results are sumjnarized in Figure 5 and numeri- 
cally in Appendix D. 



27 





N 



O 



:=) 

O' 

w 

P-. 



28 



Figure 5. Magnetic spectra computed from SOWM data applied to equations 25 and 26 
5a, data generated for 270° aircraft motion. 

5b, data generated for 000° aircraft motion. 



V. DISCUSSION 



Baker and Graefe (1968) showed experimentally that the 
theory reviewed in Section II describes the magnetic field 
induced by ocean surface gravity waves. Comparison of the 
aircraft magnetic power spectra, Figure 4, and the expected 
magnetic power spectra computed from SOWN data, Figure 5, 
further verifies the electromagnetic theory described above. 

Figure 6, which is the superpositioning of Figure 4a on 
Figure 5a and Figure 4b on.- Figure 5b, illustrates that although 
the aircraft and wave data magnetic spectra are not the same 
they are similar in amplitude and frequency. The differences 
observed in the low frequency end can, in part, be attributed 
to the assumptions necessary to the formulation of each magnetic 
power spectrum and the assumption that each data source is 
stationary in time and space . 

The FNWC data, although close to the flight path, was not 
coincident with the aircraft track (see Figure 2). Furthermore, 
determining the barometeric surface pressure field over data 
sparse areas, such as the Gulf of Alaska caused inaccuracies 
in the SOWM wave height, frequency and directional information, 
which could have contributed to some of the differences observed. 

When studying the magnetic spectrum computed from the air- 
craft data it cannot be assumed that the only source of noise 
was the ocean gravity waves. Although no significant electro- 
magnetic micropulsations in the 0 to 2 Hz frequency band, were 
expected, they cannot be ignored. 
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Figure 6. Superposition of magnetic spectra observed from aircraft ( ) and magnetic 

spectra computed from SOWM data applied to equations 25 and 26. 

6a is data from fii'’st leg of aircraft flight. 

6b is data from second leg of aircraft flight. 



Production of magnetic noise by the aircraft's motion and 
onboard equipage was considered. Ideally, most aircraft pro- 
duced noises were compensated for in the magnetometers real 
time processing scheme and the flight was straight and level 
on each leg which would further reduce noise. However, expe- 
rience shows that the noise sources described above are not 
completely removed from the spectral analysis when long data 
records are required. 

Note in Figure 6a that both the aircraft and the analyzed 
SOWN magnetic power spectra indicate a broad coincident frequency 
band. In Figure 6b it is seen that a change in the flight path 
causes the band width of coincidence to compress, shift to 
higher frequency, and two distinct maxima are observed in the ^ 

aircraft spectrum. One of which is nearly coincident with the 
SOWM magnetic spectrum. This indicates that a significant mag- 
netic noise source is external to the aircraft and ocean surface 
wave dependent . 

Most significant, is the information implied from the analysis 
of Figure 4 and Figure 6. It is seen (Figure 4) that the co- 
herence is high and the phase relationship is good bet;'?een the 
magnetometer and the altimeter in a small frequency band (.2 to 
.7 Hz). The trend observed in the phase relationship over that 
band can be attributed to signal delays imposed by the electrical 
characteristics of each instrument. 

This highly correlated frequency band falls in the same 
band in which the ocean gravity wave noise is expected. The 
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altimeter specifications indicated that it was accurate, ^ 1 
foot and has a linear response to about 1.5 Hz. Thus it can 
be concluded that the altimeter senses the sea surface flucuations . 



meter power spectra of Figure 4. The peak is approximately 
one order of magnitude greater in the data of the first flight 
leg (Figure 4a) when compared to the second (Figure 4b). This 
peak was due to an intermittent square wave caused by malfunc- 
tion of the altimeter. The square wave does not significantly 
degrade the analysis presented above . 



Note the large spectral peak below 
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VI. APPLICATION 



It was shown in section V that the altimeter can be used 
to describe the instantaneous fluctuations of the sea surface 
below the aircraft which are directly related to a portion of 
the magnetic noise. Thus, the high correlation between the 
altimeter and the magnetometer indicates that it is possible 
to remove a portion of the magnetic noise due to ocean surface 
gravity waves . 

Current Naval ASW tactics require that the magnetic field 
information be analyzed in the real time domain. Thus, noise 
cancelling must be done real time. 

Spectral analysis or convolution filtering is not practical 
due to the* signal characteristics of the target anomaly and 
the long length of records required to compute the FFT which 
describes the low frequency signals. 

The magnetic anomaly is a single transient signal imposed 
on the total magnetic field thus making spectral analysis 
impossible. Furthermore, it is possible that the aircraft 
may change direction and flight attitude before a sufficiently 
long record is collected from the altimeter and the magneto- 
meter to determine the coefficients of the Fourier Transform 
for convolution filtering. 

In order to keep noise removal in the real time domain the 
principle of Adaptive Noise Cancelling (Widrow, et.al., 1975) 
was applied to the 10 Feb. data. 
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The adaptive noise cancelling was done by a microprocessor 
developed at Naval Oceans System Center (NOSC) and already 
used in frequency analysis at frequencies higher than those 
analyzed for ultra low frequency magnetic information (0-2 Hz). 
The altimeter (reference signal) contained noise that was corre- 
lated with noise in the magnetometer (primary signal). Each 
signal was digitized and the algorithm used computed the solu- 
tion which removed the correlated noise from the primary by 
adaptively subtracting the reference signal to converge on 
minimum power in the primary signal. 

Preliminary results indicate that this method may be useful 
for real time noise cancelling in airborne magnetic data (see 
Figure 7 and Figure 8) while the magnetic anomalies not due 
to wave motion should be relatively unaffected. 

Note in Figure 7 the relative signal levels before and 
after noise cancelling. The signal appears to be cut in gain 
approximately four to one after noise cancelling. However, 
spectral analysis (Figure 8) indicates that the only power 
reduction was in the frequency band that has been indentified 
as the ocean surface gravity wave noise frequency band. 



34 




35 



L.CMd IKAL. k at3 ^ 




Figure 8, Aircraft magnetic spectrum (A) compared to "noise 
cancelled" magnetic spectrum (B) . When altimeter 
was used as noise reference. Each spectrum has 
16 degrees of freedom. Note the band limited 
noise reduction indicated by the cross hatched 
area . 
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VII. SUMMARY 



It has been shown that ocean surface gravity waves con- 
tribute significantly to the environmental background noise 
observed in airborne magnetometry . The noise is the result 
of the magnetic field associated with electrical eddy currents 
induced by ocean wave propagation through the earth's magnetic 
field. 

The APN-194 altimeter aboard the P-3C OP.ION aircraft is 
capable of describing the wave climate below the aircraft sim- 
ultaneously with total field magnetic measurements observed 
by the AN/ASQ-8KV) magnetometer, also aboard the P-3C. Thus, 
the altimeter can be used to infer the wave noise measured by 
the magnetometer . 

Analysis demonstrated that it is possible to remove a sig- 
nificant portion of the magnetic noise due to ocean gravity 
waves by using the altimeter as a noise reference and sub- 
tracting that processed reference in real time from the mag- 
netometer time series output. 

Further research with real submarine targets is underway- 
in order that the ASW community will benefit from applications 
of this work. Development of a greatly improved magnetometry 
system is possible in the near future. 
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APPENDIX A 

MAGNETOMETER AND ALTIMETER SYSTEM CHARACTERISTICS 



I. MAGNETOMETER SYSTEM PARAMETERS - 

AN/ASW-8KV) MAGNETOMETER GENERAL CHARACTERISTICS 



Performance-Magnetometer Group 
Range 

Noise level 

Figure of merit (FOM) 

MAD bandpass (selectable) 
High Pass 

Low Pass 

Gain selection 
(full scale) 

Altitude compensation 



23,000 to 75,000 gammas 
0.016 gamma ( peak-to-peak) 
0.24 gamma 



Corner frequencies 
and 0.1 Hz 
Rejection slope 36 

Corner frequencies 
and 2.0 Hz 
Rejection slope 24 

0.1, 0.2, 0.4, 1, 
and 40 gammas 

Barometric type 
Reference altitude 
Linear range, 2 0 



0.04, 0.06, 0.03 
dB/octave 
0.2, 0.4, 0.6, 
dB/octave 
, 4, 10, 20, 

self-zeroing 

feet 
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Environmental Limits 



This equipment meets the requiremenrs of military specifica- 
tion MIL-E-5400 within the limits described below. 



Temperature 

(1) Equipment Operating: -54°C to 71°C except the 
Magnetic Detector DT- 3 2 3 /ASQ-81 ( V ) or Towed Body 
TB-6 2 3 /ASQ-81 ( V ) which have +55°C maximum 
temperature . 

(2) Equipment Nonoperating: -62°C to +85°C except the 
Magnetic Detector DT-32 3/ASQ-8K V) or Tovjed Body 
TB-6 2 3/ASQ-8 1 ( V) which have +71°C maximum 
temperature . 



Sea level to 30,000 feet 
0 to 100 percent 



Altitude 

Humidity 
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Shock 


Up to 15 g with duration of 11 
+1 milliseconds 


Vibration 

Input Power Requirements 
Magnetometer Group 


Curve IV (10 g, 5 to 500 Hz) except 
Detecting Set Control C-6983/ASQ-81( V) 
or C-9 08 6 /ASO-81 ( V ) which are curve 
III (5 g, 5 to 500 Hz) . 



115/200 volts ac 3-phase power 



Voltage limits 


108 to 118 volts 


Frequency limits 


380 to 420 Hz 


Power 

Panel Lamp Power 


150 VA maximum 


Voltage limit 


0 to 28.5 volts 


Frequency 


dc or ac 


Power 

Mangetometer Pick off Points 


10 watts maximum 


Signal pick off points in 
pins C-14 and C-15. 


the P-3C were Rack C-3, TB-433, 



39 



II. ALTIMETER SYSTEM PARAMETERS - 
AN/APM-194 RADAR ALTIMETER 

Manufactured by Honeywell, Inc,, Government and Ared Products 

Divis ion 



Specifications 



ft. accuracy to 40,000 ft 
60° beam width to 3dB point 
.1 sec RC constant 
25-35 nano second pulse 
0-25 nsec gate 
Altimeter Pick off Points 

Signal pick off in the P-3C were FS 1005, T3-416, pins 
A-5 and A-6 . 
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APPENDIX B 

SPECTRAL WAVE CLIMATE 
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Table B-2.2 SOWM Data Output, Grid Pt . 54.4N 
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Table B-2.3 SOWM Data Output 



^ \ 




Figure B.-l 



FNWC Northern Hemispheric Surface Pressure, 
OOZ 10 Feb 78. 
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* \ 




Figure B.-2 



FNWC Northern Hemispheric 850rab Heights, 
OOZ 10 Feb 78. 
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APPENDIX C 

DERIVATION OF SURFACE WAVE FREQUENCY SHIFT DUE 
TO MOTION OF AIRCRAFT RELATIVE TO THE WAVE 

The wave equation due to motion of the aircraft is, 

icot-mx^ io3't-mx 

e = e 

Where co is the angular frequency of the wave , m is the 
wave number, t is the time and x is the distance from a re- 
frence point in the direction of wave propagation. The primed 
terms are the values due to motion, 

= x+(ucosij;)t 

Where u is the aircraft velocity and ip is the angle clock- 
wise from the horizontal wave motion vector to the horizontal 
aircraft motion vector. 

Then , 

iojt-m(x+(ucosii^)t) ioj't-mx 

0 — 0 

which leads to 

^icot-muCcosii;) t _ ^ico't 

then solving for co' leads to 
CO ' = co-mu cosii; 

2 

, 2it,„ 2tt,, , , gT /o 

when CO = /f, m = /L and L = ° /2tt, 

equation 26 results 

^ = f [l - 2iTfu/ ^cosifj ] 

Co 



4 7 



f 



(26) 



APPENDIX D 

SPECTRAL VALUES OF AIRCRAFT DATA AND MODELED NOISE 

I. MAGNETIC CONSTANTS AT 52N, 156N FOR FEBRUARY 1978 
Earth's total magnetic field, = 50,807 nT 
Dip angle, I = 20° 
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THE MODELED MAGNETIC NOISE PRODUCED FROM 
THE PREDICTED SOWM SPECTRA DATA LEG 1 



FREQUENCY (Hz) 



AMPLITUDE /Hz) 



.012 
.03 
.12 
.14 
.16 
.18 
. 20 
.22 
.23 
.26 
.27 
.29 
.32 
.36 
.37 
.42 
.45 
.47 
.56 
.7 
.91 
1.1 



.0002 
.0002 
.0024 
.00 37 
.0029 
.0041 
.0014 
. 0001 
.00064 
.0015 
. 0028 
. 0102 
.0039 
.0041 
.00 76 
.0128 
. 0043 
.0054 
.00 52 
. 0016 
.0018 
.00002 



Table D-l.a The Modeled Magnetic Noise Produced From 
the Predicted SOWM Data, Aircraft Course 
270°. 
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EQUENCY (Hz) 


AMPLITUDE ( 


.11 


.00087 


.13 


.0012 


.14 


.0012 


.16 


.0037 


.18 


.0015 


.21 


.0014 


.25 


.0013 


.27 


.0057 


.28 


.0088 


. 31 


.00004 


.32 


.0037 


.36 


.0069 


.37 


.0063 


.42 


.0039 


.43 


.0096 


.47 


.0034 


.53 


.0012 


.56 


.0022 


.68 


. 0002 


.7 


.00057 


.89 


.00019 


.9 


.00018 


1.13 


.000018 



nT‘ 



Table D-l.b The Modeled Magnetic Noise Produced From 
the Predicted SOWM Data, Aircarft Course 
000 ° . 
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Table D-2 FFT Numerical Output of Aircraft Data (Course 270 ). 
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LOG SPECTRAL DENSITY /Hz 




Figure D-1 Altitude Spectrum (CS 000°). 
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Figure D-2 Magnetic Spectrum (CS 000^ 
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Figure D-3 Altitude Spectrum . 
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Figure D-4 Magnetic Spectrum 
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